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ABSTRACT 

In this paper an adaptive vibration reduction control (AVRC) design is described for a 
Stirling cryocooler combined with a passive balancer. The AVRC design was based on a mass- 
spring model of the cooler and balancer, and the AVRC algorithm described in this paper was 
based on an adaptive binary search. Results are shown comparing the baseline uncontrolled 
cooler with no balancer, the cooler with the balancer, and, finally, the cooler with the balancer 
and the AVRC. The comparison shows that it may be possible to meet stringent vibration 
reduction requirements without an active balancer. 


INTRODUCTION 

The work described herein was supported by the Space Solar Power (SSP) project to 
conduct initial cryocooler characterization studies. The SSP project proposes the use of large 
fields of solar arrays in space to produce power, which would then be transmitted to Earth 
receiving stations in microwave beams. The cables in this power distribution system are to be 
cooled at superconducting temperatures. However, cryocooler vibrations that put mechanical 
stresses on the cables could cause loss of superconductivity, resulting in catastrophic failures. 

At NASA Glenn Research Center, a cryocooler testbed has been built to address cryocooler 
performance in terms of power, temperature, vibrations, mounting techniques, and heat rejection. 
Of particular importance is the cooler’s vibration performance, since a variety of cryocooler 
applications require a reduction in vibration levels. The purpose of the vibration work at NASA 
Glenn is to establish a methodology for controlling cryocooler vibrations with a reduced parts 
count. Also, the controls must be adaptive because vibration dynamics are expected to change 
over time. 
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Another important consideration for reducing transmitted vibration forces is the selection of 
cryocooler mounting. Work is also in progress on designing and investigating cryocooler soft 
mounting techniques, but this area of study is beyond the scope of this paper. The cooler used in 
this testbed (Global Cooling Manufacturing Co., Athens, Ohio) can operate only down to about 
-120 °C. However, in terms of vibrations, the dynamics of this cooler are similar to coolers 
operating at superconducting temperatures. 

This paper describes the cryocooler testbed at NASA Glenn and presents a derivation of the 
cooler mass-spring dynamics. An overview of the adaptive vibration reduction control (AVRC) 
design is also presented. For comparison, test data from the implemented passive balancer and 
AVRC control with baseline data from the uncontrolled cooler with no balancer are also 
presented. 


CRYOCOOLER TESTBED 

The cryocooler testbed diagram is shown schematically in Figure 1 and pictorially in 
Figure 2. As seen in Figure 2, the cryocooler is surrounded by insulation on the cold tip to 
maintain cooling, and the balancer-accelerometer transducer is mounted at the end of the exposed 
part of the cooler. The control functions and parameter monitoring are performed through a 
computer using LabVIEW (National Instruments Corporation, Austin, Texas) application 
software. One cycle of the desired waveform is loaded onto a General Purpose Interface Bus 
(GPIB) 488.2 to the arbitrary waveform generator. The waveform generator continuously outputs 
the loaded signal until one cycle of a new waveform is downloaded. The waveform generator 
starts outputting the new waveform after it completes the current cycle of the old waveform 
(at zero). This allows for a seamless transition while updating a synthesized waveform 
for control purposes. The output of the waveform generator is connected to a 
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FIGURE 2. Cryocooler testbed. 


differential power amplifier, which drives the cooler. This power amplifier was also designed and 
built for this project. The computer generates the cooler power drive signal, as well as the 
vibration cancellation signal, and the two signals are superimposed to finally drive the power 
amplifier. Acceleration measurements are taken with the accelerometer amplifier from the 
accelerometer transducer and used for vibration control. The passive balancer is attached to the 
cryocooler to dampen vibrations. Other parameters being monitored include temperatures, 
current, voltages, and chiller water flow. With soft mounting, the cryocooler rests on a vibration 
isolation table to isolate it from facility-induced vibrations, as shown in Figure 2. 


CRYOCOOLER MASS-SPRING MODEL 

In this study, a two-mass-spring system model for the cryocooler-balancer was developed to 
gain insight on the general system dynamic behavior for vibration control, transmitted force, and 
the design of a soft-mounting cryocooler assembly. 

The cryocooler-balancer two-mass-spring system shown in Figure 3 can be mathematically 
described with the following system of differential equations [1]: 


m l x l + cx, + k { x, - k 2 (x 2 - x, ) = F 

(1) 

m 2 x 2 + k 2 (x 2 — x, ) = 0 

(2) 


where m\, mi are the cryocooler and balancer mass; xi, xi are the cooler and balancer translation 
distance in the cooler axial direction; k\, ki are the cooler to the external interface and the 
balancer to the cooler spring constants; c is the cooler to the external interface damping; and F is 
the cooler unbalanced force (typically sinusoidal). 

This unbalanced force has its fundamental component at the power drive frequency with its 
respective higher order harmonics. Computing the Laplace transform of 
equation (2) gives 

X 1 (s) = _col_ 

X^s) s 2 +co; w 
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FIGURE 3. Cooler-balancer mass-spring system. 


where co ? = 
transform 


2 = Jk 2 / m 2 . Dividing both sides of equation (1) by k u computing its Laplace 
, and substituting the expression of X 2 ( s ) from (3) gives 


^-5 2 X,(s) + £ S*,(s) + [ l + £. k(5)- ( Y |J ? y |W = f 

k i k i *1 I s- + co; *, 


( k 2/ k \ 


(4) 


where. F 0 = ni\g. Substituting the expressions co x = ~Jk t 7m l and the deflection X 0 = F 0 / &i in 
equation (4), multiplying both sides by the expression (s 2 + (oj) , and grouping similar terms gives 


X l (s) _ 

X « (I> 


C0 2 (s 2 +( o\) 


( m \ 
co; +(Q~ +~^co; 

m l - 


(5) 


S~ + COjS + (C0 X (0^ ) 2 

m x 1 1 ■ 


X 2 0)_X,(5) X 2 (s) 


co 2 co 2 


x 0 (s) X 0 (s)X x (s) 

s + — s + 

m, 


( 0 ; +co? + ^-co 2 Is 2 +— a> 2 S + (co,co 2 ) 2 

■ m i ~) m x 


From Figure 3, the transfer function of the transmitted force F, over F 0 is 


( 6 ) 


Wj) 


H m 

F 0 &,X 0 (j) 


■5 + 1 


X x (s) 
x 0 (s) 


, which is the same as equation (5). 


(7) 


The (cs / ki + 1) = 1 at the frequencies of interest and for the damping factors anticipated. 

Since most of the interest lies in the frequency response of F, / F 0 = X\ ! X Q , parametric 
studies involving the gain of this transfer function will be performed here, varying themass ratio 
of m\ / m 2 , the damping c, and the natural frequency of the cooler mass-spring system. The 
transfer function in equation (5) consists of 4 poles and 2 zeros. The combined 
two-mass-spring system causes the two peaks (2 double poles) to separate from their respective 
one-mass-spring natural frequencies of/i = 40 Hz and/ 2 = 60.8 Hz, respectively. The double zero 
remains at the balancer resonant frequency of f 2 . From Figure 4(a), it appears that the smaller the 
ratio of cooler to balancer mass, the further apart the separation of these frequencies on ei the r 
side of the balancer resonance. The concern was that the higher frequency peak might get into the 
region of the second harmonic of the cooler drive frequency of 60 Hz, causing amplification of 
the vibrations at the second harmonic. Examination of Figure 4(a) for reasonable mass ratio 
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eliminates this concern. The cooler mass ratio in our setup is on the order of 10. Varying the 
cooler damping constant in Figure 4(b) changes the overshoot of the low-frequency peak but has 
no other effect. Lowering the natural frequency of the cooler in Figure 4(c) causes the low- 
frequency peak to move down in frequency, and it apparently also changes the damping of this 
peak. 

It is evident from Figure 4(a), (b), and (c) that the transfer function F, / F 0 has its maximum 
attenuation at the resonant frequency of the balancer fi. This means that the fundamental of the 
transmitted force will be at its lowest when the cooler is operated at the balancer natural 
frequency. The attenuation also increases when the natural frequencies of the cooler and the 
balancer are further separated. However, there are practical considerations in separating these 
frequencies such as increased balancer mass or intrusion into launch vibration frequencies. 
Cooler testing confirmed the model predictions that the lowest vibrations are realized when the 
fundamental drive frequency of the cooler is chosen at the balancer resonant frequency. In this 
setup, the force, Fo, was calculated indirectly from acceleration measurements. However, 
decreasing Fo ultimately also decreases the transmitted force, F t , since F, = k\X\ and X) increase 
as Fo increases, with k\ constant. Additionally, testing confirmed the existence of two resonant 
peaks, one on either side of the balancer resonant frequency. 


Cooler/balancer System - X1(s)/X0(s) 



Cooler/balancer System - X1(s)/X0(s) 



Cooler/balancer System - Xl(s)/X0(s) 



FIGURE 4. Frequency response of equation (7), varying (a) the mass ratio (m, / m 2 = 1, 10, 25), Ft - 40, 
c = 157, (b) the damping (c = 60, 157, 314),/, = 40, mjm 2 = 10, and (c) the cooler mass-spring natural frequency 
(fi = 40, 25, 10) c=l51,m l /m 2 = 10. 
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Cooler/balancer System - X2(s)/X1(s) 
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FIGURE 5. Frequency response of equation (3) with / = 40,/*= 60.8, c =157, and / m 2 =10. 

Another transfer function of interest is the response of ^ 2 ( 5 ) / X[(5) in equation (3) shown in 
Figure 5. Below the resonance in Figure 5, the phase changes from 0 to 180°, which means that 
*1 and X 2 move from being in phase to opposing phase in this region. However, Figure 5 shows 
that the transmitted force is attenuated the same way on both sides of the resonance. This may not 
be exactly what would be expected, as in the case of the opposing coolers, where maximum 
vibration reduction is achieved when the two pistons are moving 180° out of phase. 


ADAPTIVE VIBRATION REDUCTION CONTROLS DESIGN AND TESTING 

The vibration cancellation signal is synthesized according to the sensed vibrations and 
superimposed on the power drive signal to the cooler. No active balancer is involved in this 
process. The Lab VIEW control development software is used for signal sensing and control. The 
hardware control setup is described in the cryocooler testbed section. 

The AVRC algorithm is based on an adaptive binary search whereby the algorithm either 
keeps stepping in a certain direction to minimize a function or performance objective (as long as 
better results are obtained) or halves the step size and reverses direction 
(if worse results are obtained). The minimum step size is limited to keep the algorithm 
continuously active. Figure 6 shows the control diagram of the AVRC binary search algorithm. 
Unlike the active balancer control approach [2], the AVRC algorithm does not depend on online 
model identification, which adds complexity. Acceleration vibrations are sensed, and fast Fourier 
transforms (FFT’s) are performed to extract the harmonic spectral content that is needed to 
synthesize the vibration cancellation signal. Before the algorithm starts, the amplitude of the 
fundamental waveform is selected. The algorithm adjusts the fundamental frequency until 
vibrations are minimized at the resonant frequency of the passive balancer. Vibrations are 
minimized at the balancer resonant frequency because the transmitted force is at its minimum at 
that point, as described in the previous section. Next, the AVRC algorithm selects the amplitude 
of the second harmonic vibration cancellation signal and starts adjusting its phase until it 
minimizes vibrations for this preselected second harmonic amplitude. The AVRC then adjusts 
the amplitude of the second harmonic vibration cancellation signal, again until vibrations are 
minimized. This process is followed for control of vibrations from the third harmonic up to the 
sixth harmonic. Finally, the algorithm loops back to repeat the process. 
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FIGURE 6. Binary search algorithm control diagram. 


Figures 7, 8, and 9 show, respectively, the spectral density of the peak acceleration with no 
balancer, with the balancer driving the cooler near the balancer resonant frequency, and with the 
balancer and the AVRC. The base in these figures is 20 logio and 1 g acceleration equal to 0 dB. 
The noise floor was somewhere around -60 dB. With the balancer in Figure 8, vibrations at the 
fundamental frequency decreased by approximately 15 dB. With the AVRC, active vibrations 
improve by another approximately 15 dB or 40 times overall, considering that in Figure 9, the 
power level was 38 percent greater and vibrations were found to scale proportionally to the 
power level. More important, the controller found the optimum drive frequency for operation. 
Without adaptive control, small variations of the balancer mass-spring parameters could result in 
significant vibration levels because of the sharpness of vibration attenuation at the balancer 
resonance. The fundamental frequency in Figure 9 shows something like a double peaking. A 60- 
Hz noise internal to the accelerometer amplifier at the sensing level of approximately -30 db g 
was found to be the cause. At different phases, the noise added to the actual vibration level, 
which caused the measurement and control to jump around. The noise was separated in 
frequency from the actual vibration signal by increasing the fidelity of the FFT’s, and edited out 
using software. This produced a smooth signal at the fundamental frequency. 

It may be possible to achieve significantly greater vibration reduction by taking advantage of 
the dip in the frequency response, as shown in Figure 4. However, FFT’s take a long time to run 
(a control-sampling time of approximately 5 sec was used). Also, increasing the phase resolution 
by increasing the data-sampling rate per frequency cycle decreases the frequency resolution for a 
fixed-size FFT window. In this case, the frequency sensing resolution is approximately 1 Hz. 
Since the system has a sharp resonance at the balancer resonant frequency, a 1-Hz frequency- 
sensing resolution is not adequate. 



FIGURE 7. Spectral of peak acceleration of cooler 
(5 V, 2. 1 A, and 60.8 Hz) with no balancer. 



Frequency (Hz) 


FIGURE 8. Spectral of peak acceleration of cooler with 
balancer driving the cooler near the balancer resonant 
frequency (5 V, 2.1 A, and 60.8 Hz). 



Frequency (Hz) 


FIGURE 9. Spectral of peak acceleration of cooler-balancer 
and AVRC (5.0 V, 2.9 A, and approx 60.8 Hz). 
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( 8 ) 


The vibration force is calculated as 


F 0 = m \ a 

where a is the acceleration in g shown in Figures 7 to 9. With the cooler weight m\ ~ 2.5 kg, the 
vibration force for the AVRC case in Figure 9 at the fundamental frequency is 

F 01 = 2.5 kg(9.8m/sec 2 ) log 10 (-27/20) = 1 . 1 N (9) 

and total force generated is 

' r o=V F oi + ^+-+^s' r o, (10) 

Maximum allowable vibrations for instrument cooling applications are typically on the order 
of 0.2 N. This level of vibration was previously achieved with control of opposing cooler 
arrangements using active balancers. It may also be achievable with passive balancer-single 
cooler configurations by improving the control design and/or introducing mounting arrangements 
that further isolate the transmitted force. 

A temperature controller was also designed and integrated with the vibration controller. This 
was a proportional integral controller, with the frequency of the controller zero chosen to give a 
90-percent response of approximately 1 10 sec. This design was chosen to avoid overdriving the 
cooler for temperature control and to allow for reasonable dynamic decoupling between the 
vibration and temperature controllers. The gain of the controller was chosen to give a reasonably 
damped response. 


CONCLUSION 

In this paper, an adaptive vibration reduction algorithm was described for a cryocooler with 
a passive balancer. The results showed that vibrations were reduced approximately 
40-fold with this balancer-controller combination. Employing a vibration-sensing technique other 
than fast Fourier transforms (FFT’s) could improve this control algorithm because of the time 
required in FFT computation and the frequency resolution needed for this type of control. 
Improvements to the controls design and/or mounting arrangements that further isolate the 
transmitted force could result in vibration reduction previously obtained only with active 
balancer-opposing cooler configurations. 
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